Clock synchronization with dispersion cancellation 
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The dispersion cancellation feature of pulses which are entangled in frequency is employed to 
synchronize clocks of distant parties. The proposed protocol is insensitive to the pulse distortion 
caused by transit through a dispersive medium. Since there is cancellation to all orders, also the 
effects of slowly fluctuating dispersive media are compensated. The experimental setup can be 
realized with currently available technology, at least for a proof of principle. 
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If two distant persons want to synchronize their clocks, 
one of them (say Alice) sends a pulse, which is bounced 
back by the other (Bob). Alice measures the pulse cre- 
ation and travel time and Bob the time at which the 
bounce occurs. By exchanging their measurement re- 
sults, each of them may know the time of the other rel- 
ative to his own clock. This is the idea underlying the 
Einstein synchronization scheme 0]. Of course, if the 
medium between Alice and Bob is dispersive, then the 
pulses they exchange get distorted and the measurements 
of the bounce time at Bob's side and of the arrival time 
at Alice's side acquire an error, which must be summed 
to the intrinsic error in such kind of measurements. In 
the actual technological problem of synchronizing clocks 
of distant parties it seems that the main cause of 
errors is given by the (possibly fluctuating) dispersion 
effects of the medium through which the pulses travel. 
Franson |^ proposed a scheme capable of suppressing 
the dispersion effects to all orders, given a suitably tai- 
lored dispersive media. Steinberg, Kwiat and Chiao 
proposed and experimentally implemented a scheme ca- 
pable of suppressing to first order the effects of dispersion 
of arbitrary media. In their scheme the time-resolved co- 
herence properties of frequency entangled pulses to first 
order do not acquire any spread when traveling through a 
dispersive medium. However their scheme is not suitable 
for clock synchronization since the dispersion cancella- 
tion is present only if the time of arrival is not deter- 
mined accurately |^,^ and since only the interferometer 
arms path length difference is recovered. 

Here a modified version of the Steinberg, Kwiat and 
Chiao interferometer is presented, which allows the syn- 
chronization of the clocks of Alice and Bob without be- 
ing bothered by the pulse distortion as it travels through 
the intervening medium. This scheme is an application 
of the proposal ||^ to employ frequency-entangled pulses 
to achieve an accuracy increase in clock synchronization. 
The synchronization protocol employed is quite differ- 
ent from the Einstein clock synchronization scheme: no 
time measurement are needed, and the relative distance 
or the transit time between Alice and Bob or any disper- 
sive property of the medium play no role in the protocol. 
As in the Einstein protocol, the only (rather reasonable) 
hypothesis is that the pulse time of travel is the same 
both ways. This allows to also employ the scheme in the 



presence of fluctuations in the medium under the require- 
ment that the fluctuations have a time scale longer than 
the pulse travel time. 

The setup is realizable with currently available tech- 
nology, since it employs conventional parametric down 
converter crystals as frequency-entanglement source. 

In the first section of this paper the time synchroniza- 
tion protocol that the proposed scheme employs is de- 
scribed. In Sect. H the experimental setup and its fea- 
tures are presented. In Sect. HI the equations of motion 



of the system are solved and the hypotheses needed to 
have dispersion cancellation are given. 



I. CLOCK SYNCHRONIZATION PROTOCOL 

In this section, the protocol underlying the proposed 
experimental setup is described and compared to the Ein- 
stein clock synchronization. It is a classical protocol 
the quantum mechanical features of the setup that will be 
introduced in the following sections are employed only to 
achieve enhanced accuracy and dispersion cancellation. 

Consider the following scenario for the sake of illustrat- 
ing the method: a conveyor belt (i.e. a physical system 
in which the transit time from A to B is the same as from 
B to A) connects Alice and Bob (as in Fig. |l|). 
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FIG. 1. 



Alice pours a quantity of sand which is proportional 
to the time shown on her clock on both sides A and A' 
of the conveyor belt, i.e., on the side traveling towards 
Bob and on the side traveling towards the point M. Bob, 
at one end of the conveyor belt (point B), scoops away 
a quantity of sand proportional to twice the time shown 
on his clock. If the proportionality constant of Alice and 
Bob are the same, then the quantity of sand at point M 
will be constant in time, after an initial transient when 
they start to act on the system. Thus, it is sufficient 
to measure such a quantity in order to recover the exact 
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time difference between Alice's and Bob's clock. The pre- 
cision with which this time difference may be recovered 
depends only on the precision with which it is possible 
to measure the quantity of sand at M. 

The main advantage of this scheme over the Einstein 
clock synchronization procedure is that no time of arrival 
measurements need to be performed. In fact, the time of 
arrival measurement of a pulse has an intrinsic unavoid- 
able error dependent on its bandwidth. The conveyor 
belt synchronization strikingly allows time synchroniza- 
tion without any timing measurements. On the other 
hand, while Einstein's protocol measures the distance be- 
tween Alice and Bob (if the pulse speed is known), this 
procedure does not allow the recovery of any informa- 
tion on this distance, unless Bob stops scooping the sand 
away. 

A practical realization of this classical scheme is read- 
ily implemented. An intense continuous polarized beam 
which travels from Alice to Bob and back is employed as 
conveyor belt. Alice, on her side, rotates at positions A 
and A' the beam polarization by an amount proportional 
to the time shown on her clock. Bob rotates the polar- 
ization at position B by an amount proportional to twice 
his time, but with opposite direction. The shift in po- 
larization measured at position M allows Alice and Bob 
to determine the difference in the time of their clocks up 
to a rotation period. In the next section, a scheme is 
introduced that exploits quantum properties to enhance 
the accuracy and cancel the effects of dispersion. 



II. EXPERIMENTAL SETUP 

In this section, an experimental scheme and procedure 
to implement the protocol described in the previous sec- 
tion is described. The experimental setup is sketched 
in Fig. H It is based on the Hong, Ou, and Man- 
del (HOM) interferometer which uses the frequency- 
entangled output beam generated by a parametric down 
converter crystal. After propagating through different 
optical paths, the signal {S) and idler (/) beams are in- 
terfered at a beam splitter. By measuring the photon 
coincidence rate Pc at the output ports 1 and 2 of the 
beam splitter, one may acquire very precise information 
on the path length difference in the two arms. The preci- 
sion limit is given by the inverse of the bandwidth of the 
down-converted state (twin beams) and is, to a large ex- 
tent, independent on the precision of the photon time of 
arrival measurement. Steinberg, Kwiat, and Chiao [^|J|] 
showed that in the presence of dispersive media in one 
arm of the interferometer, it is possible to cancel to first 
order the effects of such dispersion using detectors with 
a wide integration window. Here this idea is pursued 
further. It will be shown that, in the setup proposed in 
this paper, the dispersion effects may be canceled to all 
orders even if the same dispersive medium is present in 
both arms of the interferometer. By using the scheme 



proposed here, Alice and Bob may check whether their 
clocks are synchronized, and may keep them synchro- 
nized through a feedback loop under the hypothesis that 
the drift of their clocks is sufficiently slow. 




FIG. 2. Sketch of the experimental setup for the all-orders 
dispersion cancellation synchronization. Alice on the left and 
Bob on the right are separated by a dispersive medium. The 
frequency entangled state is produced by a parametric down 
converter crystal (PDC) and is made to interfere at the 50-50 
beam splitter (BS). The coincidence rate Pc is measured as 
a function of the position of the beam splitter SI. Alice and 
Bob introduce timing information in the setup through the 
time varying delays Sl'^, 6ll, 61'^, and 5Zf . Notice that two 
delays {Sli and 51^) increase in time, while the other two ((5Zf 
and Sll) decrease. 

As shown in Fig. ^, Alice operates the parametric 
down converter and generates the twin beams she sends 
through the interferometer. The beams travel up to 
Bob's position and are reflected back to Alice. She makes 
them interfere at a 50-50 beam splitter and measures the 
coincidence rate of the photodetector clicks at the out- 
puts of the beam splitter, as a function of the optical 
path length difference SI. She obtains a flat coincidence 
rate with a very narrow di p fo r 51 — (as in Fig. ^) . 

As will be shown in Sect |l|, in order to check whether 
their clocks are synchronized, Alice and Bob must intro- 
duce time varying delays Sl^{t), Sl^^{t), Sl^{t), and Sll{t) 
both on the idler / and signal S beam. (As examples of 
varying delays consider moving mirrors or time varying 
phase shifters). These delays play the role of adding or 
subtracting the sand from the conveyor belt. In order 
to have dispersion cancellation, four delays are needed 
instead of only the two that are sufficient for the syn- 
chronization. The protocol requires these delays to be 
linked to the time measured by Alice and Bob's clocks 
respectively. Consider, for example, the case of linear 
time dependence, i.e. 

(1) 

Sl!it) = -vit-^) 5l§{t)=vit-4), 

where tg and ^re the times at which Alice and Bob 
start their clocks and v is the delay rate (e.g. the speed 
of the moving mirrors). As shown in Fig. ^, Alice places 
one of her delays at the beginning and the other one at 
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the end of the transmission line. Bob appUes both his 
delays in the middle of the two paths. We might assume 
that Alice and Bob have identical clocks and only want 
to find out the time difference r = — between them. 
Because of the introduction of the delays (|l|), the final 
path length difference as measured by the HOM interfer- 
ometer will be affected by a factor dependent on r and 
Alice can measure it by observing the shift of the dip po- 
sition in the photon coincidence rate Pc as shown in Fig. 



III. ANALYSIS OF THE EXPERIMENT 

In this section the main result is derived. It will be 
shown that the dip in the coincidence rate graph is lo- 
cated at the position SIq — Avt, which shows that by 
measuring the Mandel dip position one can recover 
the time difference between Alice and Bob's clock. There 
is no dependence on the distance between Alice and Bob 
or on any property of the intervening medium. The only 
two hypotheses that are required are a) the travel time 
of the twin beam state from Alice to Bob is the same 
as from Bob to Alice (which gives a lower time limit for 
the fluctuations of the intervening medium) and b) the 
medium acts in the same way on the beams traveling 
in both directions (which gives a lower limit for the spa- 
tial inhomogeneities of the medium: the distance between 
the beams traveling in the two directions must be smaller 
than the inhomogeneity characteristic length). 

The twin beam state at the output of the crystal (cw 
pumped at frequency 2loo) is given by the maximally fre- 
quency entangled state 



duj (f>{Lj) \uJo + w)/|wo — w)< 
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where \uj)s and \lu)i denote the signal and idler 'fre- 
quency states' {i.e. state in which there is only one pho- 
ton at frequency lo and the vacuum for all the other fre- 
quencies), and (t>{i-j) is the spectral function of the down- 
converted light, centered in w = and characterized by 
the bandwidth Auj. The coincidence rate at the photode- 
tectors is given by the Mandel formula for photodetection 
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where T is the integration time window of the detectors. 
In Eq. (^) the electromagnetic fields at time tj at the 
output of the beam splitter, assuming for the sake of 
simplicity a linear polarization, are given by 



for J 1,2 



where A is the beam cross section and Xj is the position 
of the detector. It is possible to match the output fields of 



the beam splitter with the fields before the moving mir- 
rors by applying the beam splitter transformation on the 
mirrors' reflected flelds. This latter may be obtained by 
performing some Lorentz transformations on the input 
fields: first from the source to the moving-mirror frame 
and then back to the laboratory frame. This procedure 
yields the transformations for the field annihilation oper- 
ators when the fields are bounced off the moving mirrors. 
Analyze the idler beam first. One finds that the annihi- 
lation operator a/(aj) at the crystal position is evolved 
into a'j{ui) at Bob's position (at a distance L from the 
crystal) and into a'l{uj) at the beam splitter position (at 
a distance L' from Bob), with 



a'j{uj) 



-iLulA(t^-xo/c)-L/c]+iKl{u)) 



(5) 
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]+j:K[(cj/x)+iK^(w) 



where /3 = x = jzp is associated to the Doppler shift 
introduced by the moving mirrors, and xq is the distance 
of Alice's delay from the crystal. In Eq. (^) the terms 
K( and K^j: take into account the effect of the dispersive 
medium on the idler beam on their way to and from Bob 
respectively. Notice that because of the Doppler shift in- 
troduced by the first mirror (which, for /3 > 0, reduces 
the frequency of the beam after the mirror), is eval- 
uated at Lo/x- On the other hand, because of the com- 
pensation due to the second mirror, is evaluated at a 
frequency lu. Analogous procedure applies to the signal 
beam, resulting in the operator transformations 

a'g^ui) = asiuj) X (6) 
e»'^[iT7j(*S~*!;+^)+^^^^l+««fM+«/'('^/x) ^ 

where here a'g has been evaluated after the delay 51^ . 

As shown in Fig. ^, the modes which are detected at 
the output of the beam splitter are obtained as 



ai(c^i) = ^ [laU^i] e-'"^^'/^ + a'^(^i)] 



(7) 



where SI is the delay that Alice introduces in order to 
scan the path length dependence of Pc. Replacing Eq. 
(^ and (|^) into (j|), and taking the limit T ^ oo, which 
corresponds to wide time window at the detection [^,^, 
one obtains 



PcQc duJiduj2\{0\aiiuJi)a2{uj2)\ip)\'^ 



(8) 



where the usual limit Aw <C 2ajo was employed to sim- 
plify the frequency dependence of the field (^. Assuming 
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the symmetry (f>{uj) = (j){—uj) of the spectral function, the 
matrix element in (H) is given by 



1 



ai(wi)a2(w2)|V') = 77 + ^2 - 2wo) e*"^ x 
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^2i(t^i-c^o)[i^^(to-tg)-5i/c]-iAK(t^i) 



, (9) 



where 1^ is an overall phase term that will disappear tak- 
ing the modulus, and 



IS.k{ljj) = /if (uj) — + Kj(2wo — oj) — fcf (2a;o ~ f^) 
) + K f ( 
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is the contribution of the dispersion terms. It is imme- 
diate to see that Ak vanishes altogether at all orders 
by requiring that = n^j and = k^, which can be 
satisfied by allowing the from idler beam to propagate 
at a distance less than the spatial inhomogeities of the 
medium from the to signal beam and, equivalently, by 
allowing the to idler beam to propagate near the from 
signal beam. In this case any effect of the medium will 
be erased. Replacing Eq. (^) into (||), one obtains 



duj \(t){uj)\^ 1 - cos [2- (Sh- 51) 
* c 



where 51q 



(11) 

4uT in the non-relativistic 



or 61^ 

regime. Finally, in the simple case of a Gaussian spec- 
trum |(/)(ti;)p with variance Ao;^, Eq. (Ill]) becomes 
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(12) 



which, as shown in Fig. features a dip of width Aw^^/4 
centered in SI = SIq- By measuring the dip position, one 
can recover the time difference r between Alice and Bob's 
clocks. 




FIG. 3. Plot of Pc vs. 51. The dashed line refers to the 
case V — 0, the continuous line to the case v = 50 — . In this 
graph, the Gaussian spectrum of the down-converted beams 
has a bandwidth Au) — lO^^s"^ and r = 1 ns. 

Notice that, if the parameter v is known with an ac- 
curacy Av, then, in the non relativistic regime, the error 
in the determination of r is given by 



At 



{iAujpy 



Av 



(13) 



where the first term refers to the intrinsic accuracy in 
the dip position measurement. Suppose for example that 
the clocks are initially synchronized up to r = 1 ns and 
the twin beam bandwidth is Auj 
then to achieve an accuracy At 



- IQis Hz (as in [§), 
10^^ ns one has to 



use V ~ 500 m/s with Av < 50 m/s. 

In conclusion, the proposed protocol allows Bob to 
measure the time difference r between his and Alice's 
clock, without being affected by the dispersion of the 
intermediate medium. The accuracy of the scheme is de- 
pendent only on the bandwidth Alu of the twin beam 
and on the delay rate v, namely he can recover r with an 
error ^ l/{Auj v/c). 
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